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Abstract 

We propose an approach to enhance and direct the spontaneous emission from isolated emit- 
ters embedded inside hyperbolic metamaterials into single photon beams. The approach rests 
on collective plasmonic Bloch modes of hyperbolic metamaterials which propagate in highly 
directional beams called quantum resonance cones. We propose a pumping scheme using the 
transparency window of the hyperbolic metamaterial that occurs near the topological transi- 
tion. Finally, we address the challenge of outcoupling these broadband resonance cones into 
vacuum using a dielectric bullseye grating. We give a detailed analysis of quenching and de- 
sign the metamaterial to have a huge Purcell factor in a broad bandwidth inspite of the losses 
in the metal. Our work should help motivate experiments in the development of single photon 
sources for broadband emitters such as nitrogen vacancy centers in diamond. 
OCIS codes: 160.3918, 160.4236, 270.0270 ,250.5403, 240.6680, 160.4670 



1 Introduction 

Nanoscale light-matter interactions can be tailored using plasmonic approaches to have an impact 
on quantum information processing [[T]]. It has been conclusively demonstrated that propagating 
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surface plasmon polaritons (SPPs) can possess non-classical properties such as entanglement O 
and squeezing Furthermore, inspite of the decoherence expected from losses and electron 
collisions, these properties are preserved during propagation and are manifested in the outcoupled 
photons [4J. Thus SPPs could potentially be used as a carrier of quantum information at the 
nanoscale 

One important application where plasmonics can play a key role is in the development of 
room temperature single photon sources (Tl [U. Isolated emitters like single dye molecules or 
quantum dots emit single photons however their efficiency is low for practical applications. The 
figures of merit for single photon sources are the collection efficiency which can be increased by 
directional spontaneous emission and the Purcell effect [[T6l[l7l[l8l. Deterministic out-coupling 
of single photons and tolerance to emitter positioning is another important factor affecting device 
performance lfT9l . 

At low temperatures when the linewidth of the emitter is narrow, resonant routes using micro- 
cavities can provide the Purcell enhancement. The broadband nature of many quantum emitters 
such as nitrogen vacancy (NV) centers in diamond make them unsuitable for conventional resonant 
cavity based approaches [^. The high index of diamond also leads to total internal reflection of 
photons and a poor collection efficiency (4 %) ifTOl . Therefore, the use of nanoplasmonic struc- 
tures capable of a broadband Purcell effect is necessary for efficient extraction of single photons 
([HI Hi]]. Plasmonic metamaterials which engineer the spontaneous emission can also be used to 
enhance the absorption spectrum of the isolated emitter |fT3ll. However, a significant limitation of 
any plasmonic approach is the presence of loss and non-radiative quenching [14J. Especially near 
the plasmon resonance, only a modest increase in quantum efficiency and collection efficiency is 
possible ifTSll . The key is to move away from resonant approaches and focus on low mode vol- 
ume plasmons. We emphasize that for applications requiring indistinguishable photons, photonic 
crystal and narrowband cavity approaches are ideally suited whereas plasmonic/metamaterial ap- 
proaches have an advantage for broadband applications. 
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In this paper, we show that tailored plasmonic metamaterials support collective modes which 
can efficiently channel single photons from an isolated quantum emitter to highly directional beams 
and simultaneously provide a broadband Purcell effect. Our design for the single photon source, 
compatible with emitters such as NV centers in diamond, takes into account all non-idealities in 
the structure arising from the finite patterning scale, absorption and dispersion. We also provide 
a detailed account of quenching and function away from resonance to achieve an efficient single 
photon source inspite of the losses. 



2 Quantum Resonance Cones 

Our approach relies on engineering the plasmonic Bloch modes of periodic metal-dielectric struc- 
tures. Such ID (multilayer) or 2D (nanowire) plasmonic crystals can behave as an effective meta- 
material when the lattice spacing is far below the operating wavelength [|20l[2ll[23. We consider 
a ID multilayer metal-dielectric lattice as shown in Figjl^. The short range propagating surface 
plasmon polaritons on each metal-dielectric interface couple leading to Bloch modes with unique 
properties not available in conventional photonic crystals. They are described by effective medium 
theory (EMT) in the metamaterial limit by a homogeneous medium that has a metallic dispersion 
along one direction but a dielectric dispersion in the perpendicular direction. 

The multilayer structure behaves as an effective metamaterial slab with an extremely anisotropic 
dielectric tensor e = diagfey, ey, e^], where the directions are parallel and perpendicular to the 
above mentioned layers. Extraordinary plane wave propagation in uniaxial anisotropic media is 
governed by the dispersion relation 

kl/e^ + kl/e\\ = {u/cf . (1) 
which describes an open hyperboloid (FigjlTb)) when e\\e^ < 0. These artificial media are known 
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as hyperbolic metamaterials (HMM) ||23]| . They have the property of large momentum bulk propa- 
gating waves (high- A; modes with unbounded magnitudes of kx and k^) which arise due to surface- 
plasmon-polariton Bloch waves in the plasmonic crystal. For ordinary materials such as glass with 
e\\ej_ > 0, the dispersion relation describes a bounded sphere and the magnitudes of kx and kz have 
an upper cut-off. Above this cut-off, waves are evanescent and simply decay away. 

The Poynting vector in a medium is related to the normal vector to the dispersion relation (blue 
arrows m Figlltb)). The Poynting vector for various plane waves in the HMM, given by ^ = 
kx/ {koe±)x + kz/ {e\\ko)z, where kx and k^ are related by the dispersion relation, eq. [TJ lie within 
a narrow region known as the resonance cone. Thus energy flow due to the bulk plasmonic Bloch 
modes in this metamaterial is inherently directional. The half angle of the plasmonic resonance 
cone is given by ll24l 



tan 9rc 




Furthermore, for all the waves in the medium with wavevectors along the asymptotes of the hyper- 
bola, the Poynting vectors point in the same direction. Since there are infinitely many waves (in the 
EMT limit) with wavevectors along the asymptotes of the hyperbola, a spatial crowding of Poynt- 
ing vectors is expected in preferred directions (Fig|T]^b)). This phenomenon has been observed 
in anisotropic plasmas [[24l|25l and emulated in microwave metamaterial circuits (261. Here, we 
consider the plasmonic equivalent of resonance cone behavior at optical frequencies. Below, and 
in the following sections we show that spontaneous emission from an isolated emitter is enhanced 
and directed into sub-diffraction resonance cones. The enhanced emission leads to single photon 
resonance cones which can be outcoupled to vacuum, thus allowing the study of non-classical light 
propagation in metamaterials and opening the route to quantum applications. 

We consider the case of an emitter such as a quantum dot or dye molecule with a dipole tran- 
sition, placed within the practical multilayer metal-dielectric metamaterial (Fig(T]^a)). Preferential 
emission into the resonance cone is observed for a dipole emitter placed within a 30 nm Ti02 layer 
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surrounded on top and bottom by alternating layers of silver and Ti02. The thickness of the silver 
layer and Ti02 layers are 10 nm and 30 nm respectively. There are 5 Ag layers and 4 Ti02 layers on 
each side of the embedded layer. The resonance cone emission in the finite multilayer realization 
matches quite closely with the high-/c limit predicted by effective medium theory for an ideal hy- 
perbolic metamaterial. The Poynting vector is non-zero only outside the resonance cone {9 > 9rc) 
and energy tends to flow along the surface of the resonance cone 6 ^ 6rc (Fig|TJc)inset) 



3 Decay Rate Enhancement 

Along with the directional nature of radiation, we show that a large Purcell factor in a broadband 
range is possible especially when the emitter is embedded inside the metal-dielectric multilayer 
structure. This enhancement is due to the combination of three effects: a) strong overlap of the 
plasmonic Bloch modes with the emitter b) presence of a large number of such modes and c) sub- 
wavelength confinement of the emitter. The plasmonic Bloch modes of the multilayer structure 
have a higher momentum than conventional SPPs even away from resonance. This leads to broad- 
band enhancement in the local density of states and these modes are a new radiative decay channel 
for embedded dye molecules or quantum dots [|27l|28l[29l|30l. This initial research into hyperbolic 
metamaterials has shown their potential for radiative decay engineering and here we consider the 
transmission of these modes for applications such as single photon sources. 

In the quantum mechanical approach, the decay rate of a single emitter is given by 

• n) (3) 

where (5(ro, ro, uj) is the Green's tensor of the electromagnetic wave-equation llSTIl . evaluated at the 
location of the quantum emitter with dipole moment jl. Here we use the semi-classical approach 
developed by Ford and Weber fT45 which treats the quantum emitter as a radiating point dipole. 



jl • Im 



G(ro,ro,cj) 
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Figure 1: (a) A multilayer stack of Ag/Ti02 with 10/30 nm layer thicknesses behaves like a homo- 
geneous metamaterial slab with hyperbolic dispersion (ey < 0, > 0) above A ^ 720 nm. The 
emitter is embedded symmetrically in a 30 nm layer of Ti02. (b) The direction of power flow is 
normal to the isofrequency surface. In HMMs, the power flow of all high-/c states tends to bunch 
and point in the same direction, thus forming resonance cones, (c) A quantum dot embedded in 
a HMM emits preferentially into high-/c states within the resonance cone (A = 800 nm). The 
Power Density and Poynting Vector are are normalized by the total time averaged power emitted 
from the oscillating point dipole. The grey arrows show the asymptotic direction of power flow for 
extremely high-/c states for an ideal HMM. 
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In the weak coupling limit, the quantum mechanical approach and the semi-classical theory are 
equivalent. 

Using this approach, we find for a point dipole with dipole moment jl placed a distance d above 
(or embedded in) a HMM planar slab, the decay rate enhancement (relative to vacuum) to be 



/3 = r/ro = (l-r?)+ryRe 



p(A, rf, k)dk\\ 







(4) 



where k\\ is the wavevector parallel to the interface of the HMM, and rj is the intrinsic quantum 
yield. Here p(A, rf, k) is the wavevector resolved local density of states (WLDOS) seen by the 
emitter, normalized by the LDOS of vacuum. The WLDOS takes into account the various modes 
of the system which is captured by the angular reflection spectrum. 



^/if [(1 + .•<•>)*:? - (1 - rC'jfcJ] + ij^il + rW)A|| (5) 



where r^^^ and r^^^ are the reflection coefficients for s- and p-polarized light respectively, ki is the 
magnitude of the wavevector in the medium where the dipole resides ki = y/e{uj/c). Finally k^ is 



determined from the dispersion relation where the dipole resides kz = ^k\ — fcy . 

Figure [2]^a) shows the large total decay rate enhancement ;5 = F/Fo of an emitter embedded in 
a Ag/Ti02 multilayer slab across the region of hyperbolic dispersion (A > 720 nm) (Fq is the decay 
rate in a homogeneous slab of Ti02). The Green's tensor calculation of (5 shows excellent with 
the full wave FDTD calculation. The monotonically increasing behaviour of /5(A) is explained by 
noting that the coupling strength between the emitter and the HMM modes increases at shorter 
interaction distances, d/\. In figure |2]^a) d is fixed, however as A increases the interaction distance 
decreases; furthermore, material dispersion plays a strong role in the exact nature of the /3(A) 
scaling. 

This enhanced spontaneous emission is channeled to multiple modes all of which carry energy 
along the above mentioned resonance cones. We therefore analyze the modal distribution of the 
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Figure 2: (a) A large decay rate enhancement is predicted across the region of hyperbohc dis- 
persion. (b)The W-LDOS (normahzed by vacuum) available to a quantum emitter embedded in 
a multilayer realization of a HMM show that emission occurs into bulk waveguide modes of the 
HMM slab in the region of hyperbolic dispersion, (c) The transmission coefficient (computed us- 
ing a transfer matrix method for a plane wave launched from the embedded layer) at the top of the 
multilayer slab shows that the HMM waveguide modes propagate to the top of the HMM slab and 
can be outcoupled. 
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spontaneously emitted radiation using the wavevector resolved local density of states (W-LDOS) 
|[32l . This separates the density of states available to the emitter within the multilayer structure 
according to the modes characterized by the wavevector parallel to the layers (kx ). 

Figure ^h) shows the enhanced W-LDOS available to a single emitter embedded in a 30 nm 
thick layer of Ti02 clad by realistic finite slabs of Ag/Ti02 multilayer HMM (see figure [if a)). 
We identify the series of bright vertical bands as high-/c HMM waveguide modes that arise from 
coupled surface-plasmon-polaritons, and identify the last bright line as the gap mode. When the 
embedded layer thickness is kept above ^ 30 nm, emission into the high-A: modes of the HMM 
and hence resonance cones dominates the total decay enhancement /3. When the embedded layer 
is smaller, emission into the gap plasmon mode and emission quenching begins to dominate the 
total decay enhancement /3. 

Two severe limitations on plasmonic approaches in general are the finite absorption length in 
the structure and non-radiative quenching. To study the propagation length of high-/c modes in the 
hyperbolic metamaterial we plot the transfer function of a practical structure. Figure [2f c) shows 
the transmission coefficient at the top of the multilayer HMM slab, and shows that the high-/c 
HMM waveguide modes propagate a distance of L 170 nm to the edge of the slab, and therefore 
can be out coupled. A detailed account of non-radiative decay and quenching is provided later 
in this paper. We note that functioning away from resonance and optimizing the embedded layer 
thickness can lead to a large fraction of spontaneous emission directed into the high-/c modes and 
hence resonance cones. 

We emphasize that single photon sources based on the HMM high-/c modes are fundamentally 
different than those based on MIM plasmonic modes slow light modes lf35]| . localized 

plasmons [[36l[37]| and slot waveguide modes [(38]| . In MIM and slot waveguides, energy flow is 
along the MIM interfaces; while in HMM systems energy flows through the metal-dielectric multi- 
layers at oblique angles. Furthermore, by examining the band structure in the W-LDOS calculation 
we observe that the group velocity Vg = duo/ dk does not vanish and the HMM waveguide modes 
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are therefore, not slow light modes. 



4 Near field scafing laws 

Peaks in the WLDOS indicate dominant decay channels of the emitter, and as shown below, when 
the emitter is in the near-field of a HMM slab, it is dominated by unique high-/c propagating 
waveguide modes. In the following sections, the interaction distance d dependence of /3 and p for 
an effective metamaterial and a physical realizable structure are investigated. To help elucidate the 
dependence on the interaction distance we compute the power law dependence of the decay rate 
enhancement on the distance n = 

4.1 (3{d) for Effective HMM Slabs 

In this section, we treat the HMMs using effective medium theory. The systems studied are shown 
schematically in figure [3]^a,b). To compute the reflection coefficients required to calculate /3 and p 
for each of these two geometries, the fresnel reflection coefficients for plane waves incident on a 
finite slab are used. For the case of a dipole embedded between two finite HMM slabs the reflection 
coefficients are computed using plane wave reflection coefficients in the tunnel junction geometry 
which can found in lfT4ll . 

The calculated /3 and n are shown in figure [3]^c,d) for a vertically oriented dipole at a wave- 
length of A = 900 nm and for 370 nm thick metamaterial slabs consisting of Ag/Ti02 with a 
metal filling fraction of 0.25. Realistic loss and dispersion were used for the effective medium 
computation [|42l|43|]. The metamaterial is predicted to be type 2 (ey < 0, > 0) above A ^ 720 
nm; the predicted permittivity dispersion of this structure are shown in figure [8^. At A = 900 
nm, the dielectric response of the parallel and perpedicular compenent of the dielectric tensor are 
e|| = -4.3 + zO.4 and = 11.0 + i0.07 

Due to confinement effects, the embedded dipole has a larger enhanced decay rate than the 
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dipole simply above the effective medium HMM slab. Furthermore, in the limit of the dipole lo- 
cated very close to the HMM interface(s) d « A/y^, the decay rate enhancement is proportional 
to the inverse cube of the interaction distance /3 ^ d~^. This extreme near-field behaviour is at- 
tributed to emission into the unique high-/c modes of the HMM and emission into lossy surface 
waves (quenching) which is discussed in more detail in the following sections. For large interac- 
tion distances where d ^ A/^e, plane wave interference/field intensity effects govern the decay 
rate enhancement; this is discussed in more detail in the following sections. 

4.2 (3{d) for Practical Realization of HMMs: Multilayer Slabs 

In this section, we treat the HMMs using a physical realizable multilayer structure. The systems 
studied are shown schematically in figure |4]^a,b). To compute the reflection coefficients required to 
calculate /3 and p, a transfer matrix method is used. The transfer matrix method can be generalized 
in the tunnel junction geometry, where the dipole is embedded between two multilayer stacks. 

The calculated /3 and n are shown in figure Uc,d) at a wavelength of A = 900 nm and a 4.5 
period Ag/Ti02 multilayer with thicknesses 10/30 nm respectively. Realistic loss and dispersion 
were used for the calculations (421 SSJ. At A = 900 nm, Ti02 is lossless and the loss in the 
structure is due to silver which has CAg ^ — 40 + il.7 The metamaterial is predicted to be type 2 
(e|| < 0, e± > 0) at this wavelength; the predicted effective medium dispersion of this structure 
can be found in figure [8]^a). 

We observe that in general, the embedded dipole has a larger enhanced decay rate than the 
dipole simply above the multilayer. For large interaction distances where d ^ A/y^, plane wave 
interference/field intensity effects govern the decay rate enhancement; this is observed in the low-A: 
LDOS shown in following sections. When the dipole located extremely close to the HMM interface 
d « A/v^, the behaviour of the dipole located above the multilayer and the dipole embedded 
between the multilayers is different: there are different modes available for dipole to decay into in 
each geometry. This is elucidated through the high-/c LDOS discussed in the following sections. 
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5 Quenching 

The wave vector resolved density of states p(A, rf, k) provides insight into the dominate decay 
channels of the emitter and helps explain the dependence of /3 on the interaction distance. Here 
we show the preferred decay channels of a dipole at various distances from an effective medium 
HMM and a physically realizable multilayer HMM. 

5.1 Far-Field Interference 

Figure [5] shows the calculated low-k LDOS for the four geometries considered in the previous two 
sections. We see distinct peaks in the LDOS which can attributed to propagating wave interference 
effects. When the waves interfere constructively at the location of the emitter, there is an enhance- 
ment in the electric field intensity, and thus an enhancement in the decay rate of the emitter. In the 
embedded case, the propagating wave interference is due in Fabry-Perot modes. Purely dielectric 
microresonators used to enhance spontaneous emission utilize these low-fc modes. 

As the dipole is brought closer to the metal layers (embedded layer shrunk), the coupling to the 
high-/c modes increases. For interaction distances smaller than approximately half the wavelength 
in the medium d < 1/2 X/ ^/e ^ 150 nm near-field interactions begin to dominate, and emission 
into propagating waves (low-fc modes) no longer takes place. 

5.2 Near-Field Modes: EMT 

When the dipole is located in the near field of the metamaterial, emission is directed into the high-/c 
modes of HMM slab. In the effective medium limit, there is no upper cut-off on the allowed wave 
vectors, and as the dipole is brought closer to (embedded layer shrunk), extremely larger wave 
vector modes become available to the emitter. At very small distances, d/{X/y/e) <C 1 emission 
into these infinite amount of high- A; propagating modes of the HMM contribute to the enhanced 
total decay rate as the inverse cube of the interaction distance /3 ^ 1/d^ [39]. However, there 
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is a competing decay channel, namely quenching. If the emitter is brought extremely close to 
the interface then emission into lossy surface waves begins to dominate the emission into high-/c 
HMM modes, thus /5 is mainly due to non-radiative decay. This is demonstrated in figure [6] which 
shows the WLDOS p(A, rf, k) for a perpendicular oriented dipole placed above a HMM slab and 
embedded between two HMM slabs. The geometry is the same as that used in figure [3] As the 
emitter is brought closer to the HMM slab (embedded layer shrunk) coupling occurs to increasingly 
larger high-/c propagating modes. These modes are interpreted as confined waveguide modes of 
the HMM slab. In the LDOS, these modes are manifested as distinct, relatively narrow peaks. 
When the emitter is at extreme near-field distances quenching begins to dominate the total decay 
enhancement /3. Quenching (emission into lossy surface waves) appears as a broad peak in the 
LDOS. 

By separating the LDOS integral which determines the decay rate enhancement (equation |4]) 
into the two dominant decay channels, the decay rate enhancement due to the HMM waveguide 
modes (HMM) and the lossy surface waves (LSW) can be estimated from /3 ^ Phmm + f^LSw- 
Such a separation is possible by calculating the area under the WLDOS corresponding only to 
HMM modes or to LSW. By examining figure [6] we conclude that at optical wavelengths and at 
intermediate near-field distance d^5 — 20 nm, emission into high-/c HMM waveguide modes and 
hence resonance cones dominates in the EMT limit. For interaction distances ofd<5 nm, the 
emission is quenched, and the decay enhancement is largely unusable. 

5.3 Near-Field Modes: Multilayer Realization 

In the effective medium limit there is no upper cut-off on the transverse wave vector allowed by 
the medium. In a physical realizable system, the maximum transverse wave vector k\\ that can be 
supported is inversely proportional to the unit cell size 1/a. As a result, there is a finite number of 
high-/c modes supported by the multilayer HMM structure. This is demonstrated in figure[7]which 
shows the LDOS dependence on the interaction distances d for the two geometries shown in figure 
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As the emitter is brought closer the multilayer slab (embedded layer shrunk), additional high- 
k modes of the HMM appear to the emitter. We observe the presence of the finite period cut- 
off, above which no high-/c HMM waveguide modes are present. The onset of quenching (broad 
smooth peak in LDOS) occurs at much closer distances in the physical multilayer system than in 
the EMT approximation (compare figures [6] and [7]). In the case of an emitter embedded between 
two multilayer HMMs, we see the presence of a gap plasmon mode (metal-insulator-metal plasmon 
mode) when the embedded layer thickness (gap) is smaller than about 15 nm. The gap mode can 
identified by a large sharp peak in the LDOS. 

Similar to above, the LDOS integral which determines the decay rate enhancement f3 (equation 
|4]) is separated into the three dominant decay channels: high-A; bloch plasmon modes (HMM), the 
gap plasmon mode (GPM), and lossy surface waves (LSW). The decay rate enhancement from 
each can then be determined from /3 ^ I^hmm + I^gpm + I^lsw Utilizing figure [7| we conclude 
that at optical wavelengths and at intermediate near-field interaction distances of d ^ 10 — 20 
nm, emission into high-k modes {(^hmm) is much greater than the other two contributions. These 
high-/c modes contribute to emission into highly directive resonance cones. By moving the emitter 
closer to the multilayer HMM, emission into the gap mode and quenching begins to dominate the 
total decay enhancement. However, in contrast to the EMT case, when the emitter is embedded in 
an extremely small layer, emission into the gap plasmon mode dominates, not quenching, the other 
decay channels. This result is in agreement with [33.1 . The gap plasmon mode does not contribute 
to emission into resonance cones. 

6 Topological transition assisted pumping scheme 

We now address the critical challenge of pumping embedded emitters in metamaterials. Quite often 
the large impedance mismatch with vacuum leads to inefficient coupling of radiation to emitters 



14 



placed inside such structures. We design our metamaterial to have a transparency window at the 
pumping wavelength and a large density of states (plasmonic response) only in the wavelength of 
emission. 

Figure [8]^a) shows the effective medium dielectric response of the Ag/Ti02 multilayer. Below 
XoTT ^ 720 nm the multilayer slab is effectively a dielectric (ey > 0, > 0) and transparent to 
normal incident light. Above Xqtt ^ 720 nm the multilayer slab is a type 2 HMM (ey > 0, < 
0), where enormous increase in the LDOS yields large decay rate enhancement /3 into high-/c modes 
of the HMM. This significantly different behaviour is attributed to an optical topological transition 
(OTT) in the iso-frequency (fc-space) surface from a closed ellipsoid to an open hyperboloid at 
XoTT ^ 720 nm which drastically changes the decay channels available to the emitter (see Eqjl]) 
|[39l . The difference in the electromagnetic responses above and below the topological transition 
allows for a novel realistic pumping scheme: a single emitter can be pumped where the HMM 
is transparent, and in turn have its lower energy emission enhanced and extracted in the type 2 
hyperbolic dispersion region. 

To understand the nature of spontaneously emitted radiation we consider the case of a high 
index superstrate (e >> 1) and define the Purcell factor using only propagating waves that reach 
the far-field. As expected, without the superstrate all the modes are confined to the structure and 
do not out couple to the far-field. However a very interesting transitional behavior is seen in the 
far-field with the high index superstrate. Figure [8|3 shows the predicted far-field Purcell factor Fp 
of the emitter inside a practical metal-dielectric multilayer with a high-index superstrate (e ^ 30). 
Enhanced spontaneous emission reaches the far-field of the superstrate and there is a large increase 
in the Purcell factor exactly at the transition wavelength predicted by effective medium theory 
(Xott)' Note the broad bandwidth in which this Purcell factor occurs. 
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7 Outcoupling 

We now address the outcoupling of high-/c modes which has been the significant hmiting factor 
for all applications of hyperbolic metamaterials [[39l|28]]. Given the cylindrical symmetry of the 
dipole/multilayer system, we propose a BuUseye grating structure BOl to out couple the resonance 
cone to vacuum. Figure [9f a) shows the proposed buUseye structure etched into an 85 nm Ti02 
layer on top of the Ag/Ti02 HMM structure. Figure [9]^b) shows that the far-field Purcell factor can 
exceed 6, and the spectral location of the maximum Fp can be tuned by varying the BuUseye grating 
period A. It is worth emphasizing that the colour selective out coupling by the grating is a resonant 
phenomenon ll4T]| . however the magnitude of the Fp is due to the decay rate enhancement into the 
high-/c HMM modes (resonance cones), which is fundamentally non-resonant and broadband. 

Since the quantum emitter couples into specific metamaterial modes with a unique lateral wave 
vector kx (figure[2^), the BuUseye grating scatters these modes into well defined free space modes. 
The directivity D gives a measure of the enhanced directivity of a given emission pattern relative to 
an isotropic point source; it is defined as D{9) = P{9) / {Prad/^^) where P{9) is the power density 
at a given far-field angle 6 and Prad/^^^ is the total power radiated into the far-field averaged over 
all solid angles. Figure [9fc) shows the directivity on a linear scale for a grating period of A = 400 
nm at the wavelength of maximum far-field Purcell factor A 900 nm. We see that the BuUseye 
scatters the resonance cone into a highly directive well defined spatial mode: a thin conical shell 
with a half angle ofO ^ 34^. 

In this investigation, we have considered a single period BuUseye grating which can be used to 
outcouple a single wavelength despite the fact that the HMM structure provides a spectrally broad 
decay rate enhancement. Recent work on broadband plasmonic scattering [|44l |45]| has demon- 
strated that efficient outcoupling of plasmons across the optical spectrum can be achieved with 
chirped gratings. Modification of the BuUseye grating design can lead to broadband outcoupling 
of the single photons. The HMM approach presented in this paper can be far more advantageous 
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than microcavities for coupling with broadband emitters such as NV centers in diamond. 

8 Conclusion 

In conclusion, we have proposed a device using plasmonic metamaterials to efficiently extract sin- 
gle photons from broadband quantum emitters such as NV centers in diamond. Detailed analysis 
of quenching and non-radiative decay shows that the device can achieve a broadband purcell effect 
and directional far-field spontaneous emission despite the metal losses. Low mode volume con- 
finement of the emitter and the ability to function away from plasmonic resonances holds the key 
to efficient outcoupling of light. The proposed structure will lead to a platform for understanding 
propagation of non-classical light in plasmonic metamaterials and can be integrated with other 
plasmonic structures for future nanoscale quantum information applications. 

Acknowledgements 

Z.Jacob acknowledges E. E. Narimanov and V. M. Shalaev for discussions. This work was sup- 
ported by the Natural Science and Engineering Research Council of Canada (NSERC), Alberta 
Innovates Technology Futures (AITF), Alberta Nanobridge and Canadian School of Energy and 
Environment. 

References 

[1] Michael A. Nielsen and Isaac L. Chuang. Quantum Computation and Quantum Information. 
Cambridge University Press, December 2010. 

[2] E. Altewischer, MP Van Exter, and JP Woerdman. Plasmon-assisted transmission of entan- 
gled photons. Nature, 418(6895):304-306, 2002. 



17 



[3] A. Huck, S. Smolka, P. Lodahl, A.S. SUOOFSrensen, A. Boltasseva, J. Janousek, and U.L. 
Andersen. Demonstration of quadrature-squeezed surface plasmons in a gold waveguide. 
Physical review letters, 102(24), 2009. 

[4] Giuliana Di Martino, Yannick Sonnefraud, Stephane Kena-Cohen, Mark Tame, Aahin K. 
Azdemir, M. S. Kim, and Stefan A. Maier. Quantum statistics of surface plasmon polaritons 
in metallic stripe waveguides. Nano Letters, 12(5):2504-2508, May 2012. 

[5] Z. Jacob and V. M. Shalaev. Plasmonics goes quantum. Science, 334(6055):463-464, 2011. 

[6] Z. Jacob Quantum Plasmonics MRS Bulletin, 37(8),761-767,2012 

[7] B. Lounis and M. Orrit. Single-photon sources. Reports on Progress in Physics, 68(5): 1129, 
2005. 

[8] P. Grangier, B. Sanders, and J. Vuckovic. Focus on single photons on demand. New Journal 
of Physics, 6(1), 2004. 

[9] I. Aharonovich, A.D. Greentree, and S. Prawer. Diamond photonics. Nature Photonics, 
5(7):397-405,2011. 

[10] T. M. Babinec, B. J. M. Hausmann, M. Khan, Y. Zhang, J. R. Maze, P. R. Hemmer, and 
M. Loncar. A diamond nanowire single-photon source. Nature nanotechnology, 5(3): 195- 
199, 2010. 

[11] J.T. Choy, B.J.M. Hausmann, T.M. Babinec, I. Bulu, M. Khan, R Maletinsky, A. Yacoby, 
and M. Loncar. Enhanced single-photon emission from a diamond-silver aperture. Nature 
Photonics, 5(12):738-743, 2011. 

[12] A. Huck, S. Kumar, A. Shakoor, and U.L. Andersen. Controlled coupling of a single nitrogen- 
vacancy center to a silver nanowire. Physical Review Letters, 106(9): 96801, 201 1. 



18 



[13] K. Tanaka, E. Plum, J. Y. Ou, T. Uchino, and N. 1. Zheludev. Multifold enhancement of quan- 
tum dot luminescence in plasmonic metamaterials. Physical review letters, 105 (22): 227403, 
2010. 

[14] GW Ford and WH Weber. Electromagnetic interactions of molecules with metal surfaces. 
Physics Reports, 113(4): 195-287, 1984. 

[15] G. Sun, J.B. Khurgin, and R.A. Soref. Practicable enhancement of spontaneous emission 
using surface plasmons. Applied physics letters, 90(1 1):1 11107-11 1107, 2007. 

[16] R. Esteban, T. V. Teperik, and J. J. Greffet. Optical patch antennas for single photon emission 
using surface plasmon resonances. Physical review letters, 104(2):26802, 2010. 

[17] I. Friedler, C. Sauvan, J. P. Hugonin, P. Lalanne, J. Claudon, and J. M. Grard. Solid-state 
single photon sources: the nanowire antenna. Optics express, 17(4):2095-2110, 2009. 

[18] K. G. Lee, X. W. Chen, H. Eghhdi, R Kukura, R. Lettow, A. Renn, V. Sandoghdar, and 
S. Gtzinger. A planar dielectric antenna for directional single-photon emission and near-unity 
collection efficiency. Nature Photonics, 5(3): 166-169, 2011. 

[19] D. Englund, B. Shields, K. Rivoire, F. Hatami, J. Vuckovic, H. Park, and MD Lukin. De- 
terministic coupling of a single nitrogen vacancy center to a photonic crystal cavity. Nano 
letters, 10(10):3922, 2010. 

[20] Y. Xiong, Z. Liu, C. Sun, and X. Zhang. Two-dimensional imaging by far-field superlens at 
visible wavelengths. Nano letters, 7(ll):3360-3365, 2007. 

[21] AV Kabashin, P. Evans, S. Pastkovsky, W. Hendren, GA Wurtz, R. Atkinson, R. Pollard, 
VA Podolskiy, and AV Zayats. Plasmonic nanorod metamaterials for biosensing. Nature 
materials, 8(11):867-871, 2009. 



19 



[22] J. Elser, V.A. Podolskiy, 1. Salakhutdinov, and 1. Avrutsky. Nonlocal effects in effective- 
medium response of nanolayered metamaterials. Applied physics letters, 90(19): 191 109- 
191109, 2007. 

[23] D. R. Smith, D. Schurig, J. J. Mock, P. Kolinko, and P. Rye. Partial focusing of radiation by 
a slab of indefinite media. Applied physics letters, 84(13):2244-2246, 2004. 

[24] R. K. Fisher and R. W. Gould. Resonance cones in the field pattern of a short antenna in an 
anisotropic plasma. Physical Review Letters, 22(21): 1093-1095, 1969. 

[25] L.B. Felsen, N. Marcuvitz, et al. Radiation and scattering of waves. IEEE press Piscataway, 
NJ, 1994. 

[26] K. G. Balmain, A. A. E. Luttgen, and P. C. Kremer. Resonance cone formation, reflection, 
refraction, and focusing in a planar anisotropic metamaterial. Antennas and Wireless Propa- 
gation Letters, IEEE, 1(1): 146-149, 2002. 

[27] Zubin Jacob, Igor I Smolyaninov, and Evgenii E Narimanov. Broadband purcell effect: Ra- 
diative decay engineering with metamaterials. Applied Physics Letters, 100(18): 181 105- 
181105-4, May 2012. 

[28] M. A. Noginov, H. Li, Y. A. Bamakov, D. Dryden, G. Nataraj, G. Zhu, C. E. Bonner, 
M. Mayy, Z. Jacob, and E. E. Narimanov. Controlling spontaneous emission with meta- 
materials. Optics letters, 35(11): 1863-1865, 2010. 

[29] Z. Jacob, J.Y. Kim, GV Naik, A. Boltasseva, EE Narimanov, and VM Shalaev. Engineer- 
ing photonic density of states using metamaterials. Applied Physics B: Lasers and Optics, 
100(1):215-218, 2010. 

[30] I. lorsh, A. Poddubny, A. Orlov, P. Belov, and Y. S. Kivshar. Spontaneous emission enhance- 
ment in metal-dielectric metamaterials. Physics Letters A, 201 1. 

20 



[31] L. Novotny and B. Hecht. Principles ofnano-optics. Cambridge Univ. Press, 2006. 

[32] C. L. Cortes, W. Newman, S. Molesky, and Z. Jacob. Quantum nanophotonics using hyper- 
bolic metamaterials. Journal of Optics, 14(6):063001, 2012. 

[33] K. J. Russell, T. L. Liu, S. Cui, and E. L. Hu. Large spontaneous emission enhancement in 
plasmonic nanocavities. Nature Photonics, 6(7):459-462, 2012. 

[34] Y. C. Jun, R. Pala, and M. L. Brongersma. Strong modification of quantum dot spontaneous 
emission via gap plasmon couphng in metal nanoshts. The Journal of Physical Chemistry C, 
114(16):7269-7273,2009. 

[35] P. Yao, C. Van Vlack, A. Reza, M. Patterson, M. M. Dignam, and S. Hughes. Ultrahigh 
purcell factors and lamb shifts in slow-light metamaterial waveguides. Physical Review B, 
80(19): 195 106, 2009. 

[36] A. F. Koenderink. Plasmon nanoparticle array waveguides for single photon and single plas- 
mon sources. Nano letters, 9(12):4228-4233, 2009. 

[37] D. E. Chang, A. S. Srensen, P. R. Hemmer, and M. D. Lukin. Quantum optics with surface 
plasmons. Physical review letters, 97(5):53002, 2006. 

[38] Q. Quan, I. Bulu, and M. Lon\vcar. Broadband waveguide QED system on a chip. Physical 
Review A, 80(1):011810, 2009. 

[39] Harish N. S Krishnamoorthy, Zubin Jacob, Evgenii Narimanov, Ilona Kretzschmar, and 
Vinod M Menon. Topological transitions in metamaterials. Science, 336(6078):205-209, 
April 2012. 

[40] H. J. Lezec, A. Degiron, E. Devaux, R. A. Linke, L. Martin-Moreno, F. J. Garcia- Vidal, and 
T. W. Ebbesen. Beaming light from a subwavelength aperture. Science, 297(5582):820-822, 
2002. 



[41] N.P. de Leon, B.J. Shields, C.L. Yu, D.E. Englund, A.V. Akimov, M.D. Lukin, and H. Park. 
Tailoring light-matter interaction with a nanoscale plasmon resonator. Physical Review Let- 
ters, 108(22):226803, 2012. 

[42] PB Johnson and RW Christy. Optical Constants of the Noble Metals. Physical Review B 6, 
1972 

[43] Xingjie Ni, Zhengtong Liu, and Alexander V. Kildishev. Nanohub PhotonicsDB: Optical 
Constants, 2010. 

[44] J.S. Bouillard, S. Vilain, W. Dickson, G.A. Wurtz, and A.V. Zayats. Broadband and broad- 
angle SPP antennas basedon plasmonic crystals with linear chirp. Nature: Scientific Report 
2:829, 2012 

[45] Qiaoqiang Gan and Filbert J. Bartoli. Surface dispersion engineering of planar plasmonic 
chirped grating for complete visible rainbow trapping. Applied Physics Letters 98, 251103, 
2011 



22 



Half-Space 

e = 9 




HMM 



e = 9 
Half-Space 



^ 10' 



10 10 
Dipole Interface Distance, d[nm] 



Half-Space 
e = 9 /it 



HMM 



e = 9 
Half-Space 




Propagating Wave 
Interference 



High-/c Modes > 
and / 
Quenching / 





10 10 
Dipole Interface Distance, d[nrr\] 



Figure 3: (a) Schematic of a dipole embedded symmetrically between two 370 nm thick HMM 
slabs composed of Ag/Ti02 with a silver filling fraction of 0.25. Type 2 hyperbolic dispersion is 
predicted above A ^ 720 nm. The coupling to the metamaterial changes as the embedded layer 
thickness is varied, (b) Schematic of a dipole placed above one of the finite HMM slabs from a. 
(c) Decay rate enhancement /3 at A = 900 nm is stronger when the dipole is embedded between 
two HMMs. (d) In the extreme near-field, the decay rate enhancement varies as the inverse cube of 
the interaction distance which is attributed to emission into high-/c HMM modes and quenching. 
In the far-field, the decay rate enhancement is governed by plane wave interference effects. 
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Figure 4: (a) Schematic of a dipole embedded symmetrically between two 4.5 period multilayer 
HMM slabs composed of Ag/Ti02 with thicknesses 10/30 nm respectively. Type 2 hyperbolic 
dispersion is predicted above A ^ 720 nm. The embedded layer thickness is varied, (b) A dipole is 
placed above one of the finite HMM slabs from a. The dipole/interface distance is varied, (c) Decay 
rate /3 calculations at A = 900 nm show that the enhancement is considerably stronger when the 
dipole is embedded between two multilayer HMMs. (d) In the extreme near-field, the decay rate 
enhancement varies as the inverse cube of the interaction distance; this is attributed to (/) emission 
into a gap plasmon mode and quenching for the embedded emitter and (//) emission quenching 
when the emitter is placed above the multilayer. In the far-field, the decay rate enhancement is 
governed by propagating wave interference effects. 
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Figure 5: Far-field LDOS: Calculated wavevector resolved LDOS at various far-field interaction 
distances at A = 900 nm show that when d > 1/2 \/ y/e 150 nm, low-/c interactions dominate. 
As the dipole is brought closer (embedded layer shrunk), emission is shifted into high-/c modes. 
The peaks in the low-A; LDOS are due to propagating wave interference. 
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Figure 6: Near-field LDOS - Effective Medium Theory: Calculated wave vector resolved LDOS 
at various near-field interaction distances at A = 900 nm show that as the interaction distance d 
is decreased, additional high-/c modes become accessible to the emitter. The HMM modes are 
recognizable as sharp peaks in p and are attributed to propagating waveguide modes of the HMM. 
When the emitter is embedded between two HMM slabs the coupling to high- A; HMM modes is 
slightly stronger and the onset of the strong near-field interactions occurs at larger d. At very 
small d, emission quenching begins to dominate. This is recognized as the wide smooth peak 
in p. Emission into the high-/c HMM modes and emission quenching cause the total decay rate 
enhancement /3 to scale as d~^ in the extreme near-field. 
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Figure 7: Near-field LDOS - Multilayer System: Calculated wave vector resolved LDOS at 
various interaction distances at A = 900 nm show that as the interaction distance d is decreased, 
additional high- A; modes become accessible to the emitter. No high-/c HMM modes exist above the 
cut-off (dashed line) related to the unit cell size. In the case of the emitter above the multilayer 
slab, the onset of emission into lossy surface waves (smooth broad peak) occurs at much smaller 
distances than in the EMT case. In the embedded case, we also see the presence of a gap plasmon 
mode, identified by the large peak which shifts with embedded layer thickness. This emission can 
be suppressed using the appropriate gap size. Again, the onset of emission into lossy surface waves 
occurs at smaller distances than in the EMT case. 
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Figure 8: (a) Parallel and perpendicular components of the dielectric permittivity predicted by 
EMT for a 10 nm and 30 nm Ag/Ti02 multilayer structure. A topological transition between 
elliptical dispersion and hyperbolic dispersion occurs just above Xqtt ^ 720 nm. (b) Using a large 
dielectric constant half-space superstrate (e >> 1), the high-/c modes of the resonance cones are 
outcoupled into well defined propagating modes in the far-field of the dielectric over a broadband 
spectral range. There is excellent agreement between the topological transition predicted by EMT 
and the large Purcell factor achieved in the practical multilayer structure. 



28 




Figure 9: (a) A sub wavelength cylindrical Bulls eye grating can be used to scatter the high-/c states 
which are intrinsically confined to the HMM slab, into vacuum propagating modes, (b) FDTD 
simulations show that the spectral location of the maximum far-field Purcell factor can be tuned 
by varying the grating period, (c) The directivity D = P{9) / {Ptot/^^) is shown at the location of 
maximum far-field enhancement (A = 400 nm; Xmax ^ 900 nm). The emission from the buUseye 
occurs into a highly directive conical shell. The small grey curve shows the directivity of a dipole 
in Ti02. Directivity is a measure of the power density directed along a particular direction relative 
to an isotropic point source emitting the same total power. 
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